Krylov subspace methods for eigenvalues with special properties and their analysis for normal matrices

نویسنده

  • Avram Sidi
چکیده

In this paper we propose a general approach by which eigenvalues with a special property of a given matrix A can be obtained. In this approach we first determine a scalar function $ : @ + @ whose modulus is maximized by the eigenvalues that have the special property. Next, we compute the generalized power iterations U,+l =if!I(A)Uj, j=O,l,..., where u. is an arbitrary initial vector. Finally, we apply known Krylov subspace methods, such as the Arnoldi and Lanczos methods, to the vector u, for some sufficiently large n. We can also apply the simultaneous iteration method to the subspace span($), . . . ,$‘} with some sufficiently large n, where x?‘) = I(I(A)xk), j = 0, 1,. , m = 1,. . . ,k. In all cases the resulting Ritz pairs are approximations to the eigenpairs of A with the special property. We provide a rather thorough convergence analysis of the approach involving all three methods as n -+ CC for the case in which A is a normal matrix. We also discuss the connections and similarities of our approach with the existing methods and approaches in the literature.

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

Implementation Aspects of Band Lanczos Algorithms for Computation of Eigenvalues of Large Sparse Symmetric Matrices

A band Lanczos algorithm for the iterative computation of eigenvalues and eigenvectors of a large sparse symmetric matrix is described and tested on numerical examples. It starts with a p dimensional subspace, and computes an orthonormal basis for the Krylov spaces of A, generated from this starting subspace, in which A is represented by a 2p + 1 band matrix, whose eigenvalues can be computed. ...

متن کامل

A Different Approach to Bounding the Minimal Residual Norm in Krylov Methods

In the context of Krylov methods for solving systems of linear equations, expressions and bounds are derived for the norm of the minimal residual, like the one produced by GMRES or MINRES. It is shown that the minimal residual norm is large as long as the Krylov basis is well-conditioned. In the context of non-normal matrices, examples are given where the minimal residual norm is a function of ...

متن کامل

Convergence Analysis of Krylov Subspace Iterations with Methods from Potential Theory

Krylov subspace iterations are among the best-known and most widely used numerical methods for solving linear systems of equations and for computing eigenvalues of large matrices. These methods are polynomial methods whose convergence behavior is related to the behavior of polynomials on the spectrum of the matrix. This leads to an extremal problem in polynomial approximation theory: how small ...

متن کامل

Numerical Solution of Linear Eigenvalue Problems

We review numerical methods for computing eigenvalues of matrices. We start by considering the computation of the dominant eigenpair of a general dense matrix using the power method, and then generalize to orthogonal iterations and the QR iteration with shifts. We also consider divide-and-conquer algorithms for tridiagonal matrices. The second part of this survey involves the computation of eig...

متن کامل

Sharp Ritz Value Estimates for Restarted Krylov Subspace Iterations

Gradient iterations for the Rayleigh quotient are elemental methods for computing the smallest eigenvalues of a pair of symmetric and positive definite matrices. A considerable convergence acceleration can be achieved by preconditioning and by computing Rayleigh-Ritz approximations from subspaces of increasing dimensions. An example of the resulting Krylov subspace eigensolvers is the generaliz...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

عنوان ژورنال:

دوره   شماره 

صفحات  -

تاریخ انتشار 2003